Application of a rank-based genetic association test to age-at-onset data from the Collaborative Study on the Genetics of Alcoholism study.

Association studies of quantitative traits have often relied on methods in which a normal distribution of the trait is assumed. However, quantitative phenotypes from complex human diseases are often censored, highly skewed, or contaminated with outlying values. We recently developed a rank-based association method that takes into account censoring and makes no distributional assumptions about the trait. In this study, we applied our new method to age-at-onset data on ALDX1 and ALDX2. Both traits are highly skewed (skewness > 1.9) and often censored. We performed a whole genome association study of age at onset of the ALDX1 trait using Illumina single-nucleotide polymorphisms. Only slightly more than 5% of markers were significant. However, we identified two regions on chromosomes 14 and 15, which each have at least four significant markers clustering together. These two regions may harbor genes that regulate age at onset of ALDX1 and ALDX2. Future fine mapping of these two regions with densely spaced markers is warranted.

permGPU: Using graphics processing units in RNA microarray association studies.

BACKGROUND: Many analyses of microarray association studies involve permutation, bootstrap resampling and cross-validation, that are ideally formulated as embarrassingly parallel computing problems. Given that these analyses are computationally intensive, scalable approaches that can take advantage of multi-core processor systems need to be developed. RESULTS: We have developed a CUDA based implementation, permGPU, that employs graphics processing units in microarray association studies. We illustrate the performance and applicability of permGPU within the context of permutation resampling for a number of test statistics. An extensive simulation study demonstrates a dramatic increase in performance when using permGPU on an NVIDIA GTX 280 card compared to an optimized C/C++ solution running on a conventional Linux server. CONCLUSIONS: permGPU is available as an open-source stand-alone application and as an extension package for the R statistical environment. It provides a dramatic increase in performance for permutation resampling analysis in the context of microarray association studies. The current version offers six test statistics for carrying out permutation resampling analyses for binary, quantitative and censored time-to-event traits.

Functional annotation signatures of disease susceptibility loci improve SNP association analysis.

BACKGROUND: Genetic association studies are conducted to discover genetic loci that contribute to an inherited trait, identify the variants behind these associations and ascertain their functional role in determining the phenotype. To date, functional annotations of the genetic variants have rarely played more than an indirect role in assessing evidence for association. Here, we demonstrate how these data can be systematically integrated into an association study's analysis plan. RESULTS: We developed a Bayesian statistical model for the prior probability of phenotype-genotype association that incorporates data from past association studies and publicly available functional annotation data regarding the susceptibility variants under study. The model takes the form of a binary regression of association status on a set of annotation variables whose coefficients were estimated through an analysis of associated SNPs in the GWAS Catalog (GC). The functional predictors examined included measures that have been demonstrated to correlate with the association status of SNPs in the GC and some whose utility in this regard is speculative: summaries of the UCSC Human Genome Browser ENCODE super-track data, dbSNP function class, sequence conservation summaries, proximity to genomic variants in the Database of Genomic Variants and known regulatory elements in the Open Regulatory Annotation database, PolyPhen-2 probabilities and RegulomeDB categories. Because we expected that only a fraction of the annotations would contribute to predicting association, we employed a penalized likelihood method to reduce the impact of non-informative predictors and evaluated the model's ability to predict GC SNPs not used to construct the model. We show that the functional data alone are predictive of a SNP's presence in the GC. Further, using data from a genome-wide study of ovarian cancer, we demonstrate that their use as prior data when testing for association is practical at the genome-wide scale and improves power to detect associations. CONCLUSIONS: We show how diverse functional annotations can be efficiently combined to create 'functional signatures' that predict the a priori odds of a variant's association to a trait and how these signatures can be integrated into a standard genome-wide-scale association analysis, resulting in improved power to detect truly associated variants.

Novel associations of UDP-glucuronosyltransferase 2B gene variants with prostate cancer risk in a multiethnic study.

BACKGROUND: We have previously shown that a functional polymorphism of the UGT2B15 gene (rs1902023) was associated with increased risk of prostate cancer (PC). Novel functional polymorphisms of the UGT2B17 and UGT2B15 genes have been recently characterized by in vitro assays but have not been evaluated in epidemiologic studies. METHODS: Fifteen functional SNPs of the UGT2B17 and UGT2B15 genes, including cis-acting UGT2B gene SNPs, were genotyped in African American and Caucasian men (233 PC cases and 342 controls). Regression models were used to analyze the association between SNPs and PC risk. RESULTS: After adjusting for race, age and BMI, we found that six UGT2B15 SNPs (rs4148269, rs3100, rs9994887, rs13112099, rs7686914 and rs7696472) were associated with an increased risk of PC in log-additive models (p < 0.05). A SNP cis-acting on UGT2B17 and UGT2B15 expression (rs17147338) was also associated with increased risk of prostate cancer (OR = 1.65, 95% CI = 1.00-2.70); while a stronger association among men with high Gleason sum was observed for SNPs rs4148269 and rs3100. CONCLUSIONS: Although small sample size limits inference, we report novel associations between UGT2B15 and UGT2B17 variants and PC risk. These associations with PC risk in men with high Gleason sum, more frequently found in African American men, support the relevance of genetic differences in the androgen metabolism pathway, which could explain, in part, the high incidence of PC among African American men. Larger studies are required.

Finding our way through phenotypes.

Despite a large and multifaceted effort to understand the vast landscape of phenotypic data, their current form inhibits productive data analysis. The lack of a community-wide, consensus-based, human- and machine-interpretable language for describing phenotypes and their genomic and environmental contexts is perhaps the most pressing scientific bottleneck to integration across many key fields in biology, including genomics, systems biology, development, medicine, evolution, ecology, and systematics. Here we survey the current phenomics landscape, including data resources and handling, and the progress that has been made to accurately capture relevant data descriptions for phenotypes. We present an example of the kind of integration across domains that computable phenotypes would enable, and we call upon the broader biology community, publishers, and relevant funding agencies to support efforts to surmount today's data barriers and facilitate analytical reproducibility.

Polymorphisms in the kinesin-like factor 1 B gene and risk of epithelial ovarian cancer in Eastern Chinese women.

The kinesin-like factor 1 B (KIF1B) gene plays an important role in the process of apoptosis and the transformation and progression of malignant cells. Genetic variations in KIF1B may contribute to risk of epithelial ovarian cancer (EOC). In this study of 1,324 EOC patients and 1,386 cancer-free female controls, we investigated associations between two potentially functional single nucleotide polymorphisms in KIF1B and EOC risk by the conditional logistic regression analysis. General linear regression model was used to evaluate the correlation between the number of variant alleles and KIF1B mRNA expression levels. We found that the rs17401966 variant AG/GG genotypes were significantly associated with a decreased risk of EOC (adjusted odds ratio (OR) = 0.81, 95 % confidence interval (CI) = 0.68-0.97), compared with the AA genotype, but no associations were observed for rs1002076. Women who carried both rs17401966 AG/GG and rs1002076 AG/AA genotypes of KIF1B had a 0.82-fold decreased risk (adjusted 95 % CI = 0.69-0.97), compared with others. Additionally, there was no evidence of possible interactions between about-mentioned co-variants. Further genotype-phenotype correlation analysis indicated that the number of rs17401966 variant G allele was significantly associated with KIF1B mRNA expression levels (P for GLM = 0.003 and 0.001 in all and Chinese subjects, respectively), with GG carriers having the lowest level of KIF1B mRNA expression. Taken together, the rs17401966 polymorphism likely regulates KIF1B mRNA expression and thus may be associated with EOC risk in Eastern Chinese women. Larger, independent studies are warranted to validate our findings.

Apoptotic variants as predictors of risk of oropharyngeal cancer recurrence after definitive radiotherapy.

Single nucleotide polymorphisms (SNPs) in the promoter region of FAS and FASLG may alter their transcriptional activity. Thus, we determined the associations between four FAS and FASLG promoter variants (FAS1377G>A, rs2234767; 670A>G, rs1800682; FASLG844T>C, rs763110 and 124A>G, rs5030772) and the risk of recurrence of squamous cell carcinoma of the oropharynx (SCCOP). We evaluated the associations between FAS and FASLG genetic variants and the risk of recurrence in a cohort of 1,008 patients. The log-rank test and multivariate Cox models were used to evaluate the associations. Compared with patients with common homozygous genotypes of FAS670 and FASLG844 polymorphisms, patients with variant genotypes had lower disease-free survival rates (log-rank p < 0.0001 and p < 0.0001, respectively) and an approximately threefold higher risk of SCCOP recurrence (HR, 3.2;95% CI, 2.2-4.6; and HR, 3.1; 95% CI, 2.2-4.4, respectively) after multivariate adjustment. Furthermore, among patients with HPV16-positive tumors, those with variant genotypes of these two polymorphisms had lower disease-free survival rates (log-rank, p < 0.0001 and p < 0.0001, respectively) and a higher recurrence risk than did patients with common homozygous genotypes (HR, 12.9; 95% CI, 3.8-43.6; and HR, 8.1; 95% CI, 3.6-18.6, respectively), whereas no significant associations were found for FAS1377 and FASLG124 polymorphisms. Our findings suggest that FAS670 and FASLG844 polymorphisms modulate the risk of recurrence of SCCOP, particularly in patients with HPV16-positive tumors. Larger studies are needed to validate these results.

Association of Common Genetic Polymorphisms with Melanoma Patient IL-12p40 Blood Levels, Risk, and Outcomes.

Recent investigation has identified association of IL-12p40 blood levels with melanoma recurrence and patient survival. No studies have investigated associations of single-nucleotide polymorphisms (SNPs) with melanoma patient IL-12p40 blood levels or their potential contributions to melanoma susceptibility or patient outcome. In the current study, 818,237 SNPs were available for 1,804 melanoma cases and 1,026 controls. IL-12p40 blood levels were assessed among 573 cases (discovery), 249 cases (case validation), and 299 controls (control validation). SNPs were evaluated for association with log[IL-12p40] levels in the discovery data set and replicated in two validation data sets, and significant SNPs were assessed for association with melanoma susceptibility and patient outcomes. The most significant SNP associated with log[IL-12p40] was in the IL-12B gene region (rs6897260, combined P=9.26 × 10(-38)); this single variant explained 13.1% of variability in log[IL-12p40]. The most significant SNP in EBF1 was rs6895454 (combined P=2.24 × 10(-9)). A marker in IL12B was associated with melanoma susceptibility (rs3213119, multivariate P=0.0499; OR=1.50, 95% CI 1.00-2.24), whereas a marker in EBF1 was associated with melanoma-specific survival in advanced-stage patients (rs10515789, multivariate P=0.02; HR=1.93, 95% CI 1.11-3.35). Both EBF1 and IL12B strongly regulate IL-12p40 blood levels, and IL-12p40 polymorphisms may contribute to melanoma susceptibility and influence patient outcome.

Impact of gene variants on sex-specific regulation of human Scavenger receptor class B type 1 (SR-BI) expression in liver and association with lipid levels in a population-based study.

BACKGROUND: Several studies have noted that genetic variants of SCARB1, a lipoprotein receptor involved in reverse cholesterol transport, are associated with serum lipid levels in a sex-dependent fashion. However, the mechanism underlying this gene by sex interaction has not been explored. METHODS: We utilized both epidemiological and molecular methods to study how estrogen and gene variants interact to influence SCARB1 expression and lipid levels. Interaction between 35 SCARB1 haplotype-tagged polymorphisms and endogenous estradiol levels was assessed in 498 postmenopausal Caucasian women from the population-based Rancho Bernardo Study. We further examined associated variants with overall and SCARB1 splice variant (SR-BI and SR-BII) expression in 91 human liver tissues using quantitative real-time PCR. RESULTS: Several variants on a haplotype block spanning intron 11 to intron 12 of SCARB1 showed significant gene by estradiol interaction affecting serum lipid levels, the strongest for rs838895 with HDL-cholesterol (p=9.2x10(-4)) and triglycerides (p=1.3x10(-3)) and the triglyceride:HDL cholesterol ratio (p=2.7x10(-4)). These same variants were associated with expression of the SR-BI isoform in a sex-specific fashion, with the strongest association found among liver tissue from 52 young women<45 years old (p=0.002). CONCLUSIONS: Estrogen and SCARB1 genotype may act synergistically to regulate expression of SCARB1 isoforms and impact serum levels of HDL cholesterol and triglycerides. This work highlights the importance of considering sex-dependent effects of gene variants on serum lipid levels.

Common genetic variation and the control of HIV-1 in humans.

To extend the understanding of host genetic determinants of HIV-1 control, we performed a genome-wide association study in a cohort of 2,554 infected Caucasian subjects. The study was powered to detect common genetic variants explaining down to 1.3% of the variability in viral load at set point. We provide overwhelming confirmation of three associations previously reported in a genome-wide study and show further independent effects of both common and rare variants in the Major Histocompatibility Complex region (MHC). We also examined the polymorphisms reported in previous candidate gene studies and fail to support a role for any variant outside of the MHC or the chemokine receptor cluster on chromosome 3. In addition, we evaluated functional variants, copy-number polymorphisms, epistatic interactions, and biological pathways. This study thus represents a comprehensive assessment of common human genetic variation in HIV-1 control in Caucasians.

A genome-wide association study identifies susceptibility variants for type 2 diabetes in Han Chinese.

To investigate the underlying mechanisms of T2D pathogenesis, we looked for diabetes susceptibility genes that increase the risk of type 2 diabetes (T2D) in a Han Chinese population. A two-stage genome-wide association (GWA) study was conducted, in which 995 patients and 894 controls were genotyped using the Illumina HumanHap550-Duo BeadChip for the first genome scan stage. This was further replicated in 1,803 patients and 1,473 controls in stage 2. We found two loci not previously associated with diabetes susceptibility in and around the genes protein tyrosine phosphatase receptor type D (PTPRD) (P = 8.54x10(-10); odds ratio [OR] = 1.57; 95% confidence interval [CI] = 1.36-1.82), and serine racemase (SRR) (P = 3.06x10(-9); OR = 1.28; 95% CI = 1.18-1.39). We also confirmed that variants in KCNQ1 were associated with T2D risk, with the strongest signal at rs2237895 (P = 9.65x10(-10); OR = 1.29, 95% CI = 1.19-1.40). By identifying two novel genetic susceptibility loci in a Han Chinese population and confirming the involvement of KCNQ1, which was previously reported to be associated with T2D in Japanese and European descent populations, our results may lead to a better understanding of differences in the molecular pathogenesis of T2D among various populations.

Genome-wide association study of Lp-PLA(2) activity and mass in the Framingham Heart Study.

Lipoprotein-associated phospholipase A(2) (Lp-PLA(2)) is an emerging risk factor and therapeutic target for cardiovascular disease. The activity and mass of this enzyme are heritable traits, but major genetic determinants have not been explored in a systematic, genome-wide fashion. We carried out a genome-wide association study of Lp-PLA(2) activity and mass in 6,668 Caucasian subjects from the population-based Framingham Heart Study. Clinical data and genotypes from the Affymetrix 550K SNP array were obtained from the open-access Framingham SHARe project. Each polymorphism that passed quality control was tested for associations with Lp-PLA(2) activity and mass using linear mixed models implemented in the R statistical package, accounting for familial correlations, and controlling for age, sex, smoking, lipid-lowering-medication use, and cohort. For Lp-PLA(2) activity, polymorphisms at four independent loci reached genome-wide significance, including the APOE/APOC1 region on chromosome 19 (p = 6 x 10(-24)); CELSR2/PSRC1 on chromosome 1 (p = 3 x 10(-15)); SCARB1 on chromosome 12 (p = 1x10(-8)) and ZNF259/BUD13 in the APOA5/APOA1 gene region on chromosome 11 (p = 4 x 10(-8)). All of these remained significant after accounting for associations with LDL cholesterol, HDL cholesterol, or triglycerides. For Lp-PLA(2) mass, 12 SNPs achieved genome-wide significance, all clustering in a region on chromosome 6p12.3 near the PLA2G7 gene. Our analyses demonstrate that genetic polymorphisms may contribute to inter-individual variation in Lp-PLA(2) activity and mass.

Assessment of LD matrix measures for the analysis of biological pathway association.

Complex diseases will have multiple functional sites, and it will be invaluable to understand the cross-locus interaction in terms of linkage disequilibrium (LD) between those sites (epistasis) in addition to the haplotype-LD effects. We investigated the statistical properties of a class of matrix-based statistics to assess this epistasis. These statistical methods include two LD contrast tests (Zaykin et al., 2006) and partial least squares regression (Wang et al., 2008). To estimate Type 1 error rates and power, we simulated multiple two-variant disease models using the SIMLA software package. SIMLA allows for the joint action of up to two disease genes in the simulated data with all possible multiplicative interaction effects between them. Our goal was to detect an interaction between multiple disease-causing variants by means of their linkage disequilibrium (LD) patterns with other markers. We measured the effects of marginal disease effect size, haplotype LD, disease prevalence and minor allele frequency have on cross-locus interaction (epistasis). In the setting of strong allele effects and strong interaction, the correlation between the two disease genes was weak (r=0.2). In a complex system with multiple correlations (both marginal and interaction), it was difficult to determine the source of a significant result. Despite these complications, the partial least squares and modified LD contrast methods maintained adequate power to detect the epistatic effects; however, for many of the analyses we often could not separate interaction from a strong marginal effect. While we did not exhaust the entire parameter space of possible models, we do provide guidance on the effects that population parameters have on cross-locus interaction.

Genetic variants in IGF-I, IGF-II, IGFBP-3, and adiponectin genes and colon cancer risk in African Americans and Whites.

PURPOSE: Evaluating genetic susceptibility may clarify effects of known environmental factors and also identify individuals at high risk. We evaluated the association of four insulin-related pathway gene polymorphisms in insulin-like growth factor-1 (IGF-I) (CA)( n ) repeat, insulin-like growth factor-2 (IGF-II) (rs680), insulin-like growth factor-binding protein-3 (IGFBP-3) (rs2854744), and adiponectin (APM1 rs1501299) with colon cancer risk, as well as relationships with circulating IGF-I, IGF-II, IGFBP-3, and C-peptide in a population-based study. METHODS: Participants were African Americans (231 cases and 306 controls) and Whites (297 cases, 530 controls). Consenting subjects provided blood specimens and lifestyle/diet information. Genotyping for all genes except IGF-I was performed by the 5'-exonuclease (Taqman) assay. The IGF-I (CA)(n) repeat was assayed by PCR and fragment analysis. Circulating proteins were measured by enzyme immunoassays. Odds ratios (ORs) and 95 % confidence intervals (CIs) were calculated by logistic regression. RESULTS: The IGF-I (CA)( 19 ) repeat was higher in White controls (50 %) than African American controls (31 %). Whites homozygous for the IGF-I (CA)(19) repeat had a nearly twofold increase in risk of colon cancer (OR = 1.77; 95 % CI = 1.15-2.73), but not African Americans (OR = 0.73, 95 % CI 0.50-1.51). We observed an inverse association between the IGF-II Apa1 A-variant and colon cancer risk (OR = 0.49, 95 % CI 0.28-0.88) in Whites only. Carrying the IGFBP-3 variant alleles was associated with lower IGFBP-3 protein levels, a difference most pronounced in Whites (p-trend <0.05). CONCLUSIONS: These results support an association between insulin pathway-related genes and elevated colon cancer risk in Whites but not in African Americans.

Distribution, ecology, life history, genetic variation, and risk of extinction of nonhuman primates from Costa Rica

We examined the association between geographic distribution, ecological traits, life history, genetic diversity, and risk of extinction in nonhuman primate species from Costa Rica. All of the current nonhuman primate species from Costa Rica are included in the study; spider monkeys (Ateles geoffroyi), howling monkeys (Alouatta palliata), capuchins (Cebus capucinus), and squirrel monkeys (Saimiri oerstedii). Geographic distribution was characterized accessing existing databases. Data on ecology and life history traits were obtained through a literature review. Genetic diversity was characterized using isozyme electrophoresis. Risk of extinction was assessed from the literature. We found that species differed in all these traits. Using these data, we conducted a Pearson correlation between risk of extinction and ecological and life history traits, and genetic variation, for widely distributed species. We found a negative association between risk of extinction and population birth and growth rates; indicating that slower reproducing species had a greater risk of extinction. We found a positive association between genetic variation and risk of extinction; i.e., species showing higher genetic variation had a greater risk of extinction. The relevance of these traits for conservation efforts is discussed.